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We introduce a method for obtaining analytic approximations to the evolution of Markovian open 
quantum systems. It is based on resumming a generalized Dyson series in a way that ensures optimal 
convergence even in the absence of a small parameter. The power of this approach is demonstrated 
by two benchmark examples: the spatial detection of a free particle and the Landau-Zener problem 
in the presence of dephasing. The derived approximations are asymptotically exact and exhibit 
errors on the per mil level over the entire parameter range. 
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Introduction. — Emergent quantum technologies em- 
ploy ever larger and noisier quantum machines, such as 
solid state devices [T], complex biomolecules [2], or high 
precision detectors [3]. While it is hard to shield such 
systems from environmental influences, physicists have 
come to think of the latter as a potentially useful re- 
source in areas ranging from quantum chemistry [IJ |S] to 
quantum optics [5HSj and quantum information [3]. 

If one seeks to exploit incoherent effects beyond the 
preparation of stationary states EHl IH] it is crucial to 
have good analytic approximations of the open system 
dynamics available, along with a fitting physical intu- 
ition. Given that many real-world applications are em- 
bedded in complex environments at high temperature, a 
Markovian treatment is here a natural point of depar- 
ture. While numerical treatments yield quantitatively 
good results [12U14) . they become generally intractable 
as a large parameter space is to be explored or the sys- 
tem size increases. Moreover, numerical approaches alone 
are generally not suited to yield deeper insights into the 
general mechanisms at work. 

In this letter we develop a framework to obtain reliable 
analytic approximations to arbitrary Markovian master 
equations, based on a formal Dyson-like expansion of the 
quantum evolution. The key ingredient is an adaptive re- 
summation of the series which optimizes its convergence 
even in the absence of a small parameter. The result- 
ing decomposition into quantum jumps and periods of 
unperturbed evolution reflects the interplay of coherent 
and incoherent dynamics in a natural way and provides a 
physically meaningful interpretation. Based on the low- 
est order terms of the optimized expansion we obtain 
highly accurate approximations for the entire dynamics. 

We start by presenting the general theory, and then 
illustrate its use and power by means of two nontrivial 
open quantum problems, reflection at a spatial detector 
[T5] and Landau-Zener tunneling with dephasing [TBI [TT] . 
The resulting approximations are asymptotically exact 
and show errors on the per mil level over the entire range 
of parameters. This substantially improves existing ap- 
proximations [ISl [inj . 



Jump expansion. — The dynamics of a Markovian open 
quantum system can be described by a master equation 
pt = C{t)pt with [20] 

N ^ 

C{t)pt = [H{t),pt]+Y.L,ptL]--[L]L,,pt] . (1) 

The operators Lj, which account for the incoherent in- 
fluence of the environment, are not unique |21| . In par- 
ticular, ([T]) is invariant under the transformation 

Lj Lj,a = Lj + otj (2) 
H{t) ^ H^{t) = Hit) - ^ ^ (a*L, - a,L]) , (3) 

with a = (ai, . . . , Oat) € (D^. 

A formal expansion of the solution pt is obtained by 
decomposing the generator ([T]) into any two parts, C{t) — 
C-o. (t) + Ja (where ol labels different possible decomposi- 
tions). Expanding the propagator Texp[/g £(t')dt'] for a 
fixed decomposition yields a Dyson-like series, the jump 
expansion [221 123j 

oo 

p, = ^p("\with (4) 

n=0 

pI"^ = / dt„z^„(t,i„)j'„pi;;-^\ (5) 

Jo 

Here Ua.{t2,ti) — Texp[J*^ £a(t')di'] propagates p from 
ti to t2, T denotes time ordering, and pj"'' = Uait, 0) Pa- 
in analogy to the usual Dyson series for unitary dy- 
namics, constituents in ([4| can be viewed as periods of 
unperturbed dynamics Ua interspersed with random per- 
turbations, or jumps, jTa- The superscript (n) labeling 
different orders in the jump expansion hence denotes the 
number of jumps associated with the respective term. 
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A natural decomposition of C is given by 

j j 
Co.{t)pt = -\ {Hf{t)p, - ptHfit)) , (7) 

where H^{t) = Ha{t) - ifi/2J2j L^^Lj^^. It has the 
special property that both the jumps i7j,a and the un- 
perturbed propagators Ua are completely positive (norm- 
decreasing) superoperators. All orders pj"' of the jump 
expansion are then valid density matrices with weights 
Wn{t) = Trpj"'' adding up to Trpt = 1. The Wn{t) may 
therefore be viewed as probabilities to register n jumps 
until t. Indeed, the jump expansion with Q, ^ pro- 
vides a physically meaningful measurement interpreta- 
tion of the open system dynamics [23 , which is essentially 
equivalent to the quantum trajectory approach [24] , 

Resummation. — In contrast to usual implementations 
of the Dyson expansion, the decomposition of C{t) into 
C-a{t) + iTct does not involve a small parameter. This 
renders the series Q with its doubtful convergence prop- 
erties of little practical use. However, we will see that a 
rapidly convergent series can be generated by choosing 
optimal decompositions a. ol C{t). 

As discussed above, the different orders in Q are char- 
acterized by weights Wn{t), which satisfy w„(0) = Sn.o- 
To maximize the contribution of the lowest order around 
t = 0, the rate of change — i9tWo(i)|t=o must therefore be 
minimal. This optimization condition determines a. at 
t = Q. For i > 0, a is a function of both time and the 
jump record = . . . ',jn,tn), which collects the 

sequence of past jumps Jj^a and their associated times. 
The weights of the n lowest orders are maximized by min- 
imizing —dt^^^QWm{t) at all times, a condition yield- 
ing the optimal Q:(i,9^"). 

Since the optimal choice a depends on the record of 
past jumps, the operator decomposition underlying the 
expansion changes from term to term. It is not at all 
obvious why such an adaptive expansion should consti- 
tute a solution of the master equation. To see that this 
is indeed the case we write the order p^"'' in terms of the 
record-conditioned branches pf^ \ 

where ' = ^a(t,fR") (i, in) J;„,a(t„,m"-i)pi!' ^- The 
time derivative of pj"'' decomposes into an and an 
"'^-conditioned term 

dtpi""^ = E ^s(t,SK")(*)Pr ^+ E '^],^Mt,'^"-^) Pt^ 

Combining the ^{"-conditioned terms in dtp^^ and in 
dtp[^^^^ yields X^fm-} ^WPt^ \ which is independent of 



a. Taking the sum over all orders, we obtain the master 
equation dtpt = ^{t)pt- 

To find the optimal adaptive decomposition Q:(i, JH"), 
we now consider the rate of change 

n 

-2Re{a*(t,^")Tr(L,pf"^)}. (8) 
The minimum is attained for the choice 

<5,(t,9^") = -Tr(L,pf")), (9) 

which ensures optimal convergence of the jump expansion 
up to n"^ oder. 

We note that this type of jump operators can also be 
found in the context of stochastic unravelings, by requir- 
ing the jumps to map into orthogonal subspaces [251 US] 
or by minimizing the entropy production of an associated 
measurement [27] . The crucial point for our purposes is 
that the adaptive decomposition defined by (|9| gives rise 
to the most efficient resummation of the jump expansion 
in the sense that even a truncated series captures the 
essential part of the entire dynamics. 

In addition to optimizing the convergence of the jump 
expansion, our goal is to make its solution analytically 
tractable. A trade-off of both aspects is obtained through 
an unbiased elimination of the state-dependence of a in 
([9]) by assuming complete ignorance, p oc 1, except im- 
mediately after a jump. This way ^ simplifies to 

o^t{3) = (10) 

The operators Li a are then piecewise constant between 
successive jumps and they depend only on the sequence 
f = Ui, ■ ■ -Jn) of previous jumps. 

Other simplifications are conceivable, 6.17. incorporat- 
ing different parts of the information (IH", t) for the adap- 
tive resummation. However, we find that ( |10[ ) proves 
surprisingly powerful in a number of applications. In 
particular, already the first two orders give rise to highly 
accurate approximations for the spatial detection of par- 
ticles and for Landau-Zener tunneling with dephasing. 

Reflection by measurement. — Our first application 
demonstrates two important aspects of the resummation: 
(i) Its rapid convergence ensures that the first terms pro- 
vide good approximations, and (ii) it offers a physical 
picture of the underlying dynamics, which can serve as 
the starting point for an analytical treatment. Specif- 
ically, we find that the resummation replaces frequent 
but insignificant projective measurements of a particle 
on the left and right half space by rare but decisive tran- 
sits across the measurement boundary. 
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FIG. 1. Reflection probability of a free particle with kinetic 
energy Ein at the boundary of a projective left-right mea- 
surement with rate 7. Op en dots: numerically exact values; 
solid line: prediction (141 of the adaptive jump expansion; 
dashed line: leading order approximation (111 [18]. The re- 



sult ([14l is asymptotically exact, with a maximal deviation 
(at Ei^hj « 0.2) of below one percent. 



Assuming that projective measurements onto the pos- 
itive and negative half-axis occur with rate 7, and that 
the outcomes are discarded, the motion of a free particle 
can be described by a master equation with H = /2m, 
Li — ^Q{-~x) and L2 = ^Jl^{x)- Equation ([Io| then 
imphes that the optimized operators ./ji,Q;(j") are iden- 
tical, alternating between ^Q(x) and y^8(— i) upon 
every jump. As a consequence, the particle must traverse 
the measurement boundary a; = between two jumps, 
with the number of transits given by the jump count. 

The particle now experiences only a finite number of 
transits before going off to the left or to the right, as can 
be confirmed by a numerical unraveling. Compared with 
the unbounded number of jumps of the original expan- 
sion this reflects the rapid convergence of the resumma- 
tion. Moreover, since every second term contributes to 
reflection (no transit, two transits, etc.), the reflection 
probability is the sum of the weights — W2n{t ~> 00) 
of all even jumps terms. 

As follows from ([7|, the lowest order contribution is 
determined by the reflection at the imaginary potential 
step —ihjQ{x). The corresponding stationary solution 

\ip^J), consisting of a plane wave to the left and an ex- 
ponential tail to the right of the step, gives the reflection 
probability 



Po(fco) = 



1 - ^1 + 2m'^i/hkl 



1 + + 2vi-ii/hkl 



(11) 



as a function of the incoming wave number fcg. The state 



after the first jump IV'!),^) is obtained by projecting 1-0^" 



(0)\ 



onto the right hand side: 



{x\ip\^^) = Q{x)M-^ exp [ix^kl + 2m-ii/h^ , (12) 



where N ~2 Im-y//cQ -I- 2mji/h. Due to the adaptive up- 
date of the jump operators, the imaginary potential step 
then changes to —01"/ Q{—x), which now acts as a scat- 
tering potential for \ip^J)- Since {k\^^^) is a Lorenzian 
wave packet, one must treat the incoming (fc < 0) and 
outgoing {k > 0) parts separately with factors Po{k) and 
1. This yields 



Pi 



dfc(e(-fc)Po(fc) + e(fc)) {k\i,l}J) (13) 



—00 
00 



+ I dt2je-'->'{^^^J\Cm):j)\il-Poiko)) 



The reflection probability of given within curly 

brackets, involves a small correction term with C'{t) — 

e{-k)Q{-x-hkt/m)e{k)+e{k)e{-x-hkt/m)e{~k), 

accounting for the norm decay during the scattering 
process at the imaginary potential. All higher order 
terms depend on the specific jump times due to the non- 
stationarity of l^^^J)- However, the reflection probabil- 



ities of IV'i"'') can be approximated by that of \ipl!'J ) m 



(13), such that we have F„+i « (1 - Pi/(1 - Po(A:o)))P« 



and hence the total reflection probability 



(l~Po(fco))' 
2 - 2Po(fco) - Pi 



(14) 



In Fig. [T] we compare this readily accessible approxi- 
mation to the results of a time-consuming simulation of 
the master equation. We find that the error remains be- 
low 1% over the entire parameter range — a considerable 
improvement of the known leading order approximation 
[TH] . It is remarkable, that this degree of accuracy is ob- 
tained with only the two lowest order reflection probabil- 
ities, highlighting the rapid convergence of the optimized 
jump expansion. 

Landau-Zener transitions with dephasing. — The sec- 
ond application is a bit more involved. It demonstrates 
that the adaptive resummation reformulates incoherent 
dynamics in a way that guides our intuition towards 
highly accurate analytic approximations. 

The Landau-Zener problem involves a pair of states 
{|1),|2)} with constant coupling a and time-dependent 
energy separation vt. Rescaling time as r = t^Jv/h one 
finds that the Schrodinger equation 



d_ 



(15) 



involves a single dimensionless quantity 5 = o? jhv, called 
adiabaticity parameter. Starting out in |1) or |2) in the 
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remote past, the probability to experience a transition is 
P((5,0) = 1 - exp(-27r(5) as r ^ oo [28]. 

The coherent evolution ( [T5| will now be supplemented 
by dephasing, as described by the operator L = y/jUz (or 
equivalently by L = 2^\2){2\). This models the influ- 
ence of a continuous energy measurement, or of rapidly 
fluctuating energy levels with dephasing rate 7. As a re- 
sult, the transition probability P(5, 7) depends on S and 
7- 

Before applying the jump expansion Q and resumma- 
tion (10), it is helpful to consider the strong dephasing 
limit. Representing p as a Bloch vector, its x- and y- 
components can be adiabatically eliminated for 7 — >■ 00, 
by setting x = y = 0. This leads to a closed evolu- 
tion equation for the z-component describing the popu- 
lations. Integration yields P{6, 00) = (1 — exp(— 47r(5))/2, 
in agreement with the Landau-Zener transfer probabil- 
ity derived in the context of rapidly oscillating energy 
levels |16) . It reflects the suppression of quantum tunnel- 
ing by frequent measurements, known from the quantum 
Zeno effect. With the exact expressions for P(i5, 0) and 
P((5, 00) at hand, we can now tackle the 7-dependence of 
the transfer probability. 

Similar to the previous application, the adaptive re- 
summation ( |10[ ) causes the jump operator to alternate 
between 2.^|2)(2| and 2y^|l)(l| upon every jump. 
Therefore, a Landau-Zener transition necessarily sepa- 
rates two consecutive jumps, which limits the number 
of jumps overall and assigns the odd jump terms to 
contribute to the transition probability, i.e. P((5, 7) — 
J2nP2n+i, with P„ = ■u;„(t ^ 00). 

The incoherent dynamics can thus be viewed as an 
alternating sequence of Landau-Zener tunnelings and 
jumps. This insight allows us to map the problem to 
a classical inhomogeneous Markov process with two al- 
ternating, time-dependent rates Xq{t) and Ai(r), corre- 
sponding to the time dependent Landau-Zener transfer 
rate and to the jump rate, respectively, 

3rP2n(T) = Ai (r)p2n-l (t) - Ao(r)p2n(T) (16) 
9rP2n+l(T) = Ao(t)p2«(t) - Ai (T)p2n+1 (t) , (17) 

with po(t) = exp(-Ao(r)), Ai(r) = X^{t)dt. 

It is natural to assume that the Landau-Zener transi- 
tions occur only during a characteristic time interval of 
length T* when the energy levels are close. Therefore, 
the probability P„ for n jumps until r = c» is equal 
to the probability for n transitions until t = r*, i.e. 

Pn = P2n-l{T*) -f P2n(T*). 

To recover the strong dephasing limit, take Ai(t) ^ 
Ao(t). After adiabatically eliminating all P2n+i an in- 
homogeneous Poisson process with rate Ao(t) is thus ob- 
tained for the P„. Comparing its probability distribution 
at time r*, P„ = Aq exp(-Ao)/n! (where A; = A,,(t*)), 
with the previously derived P((5, cx)) implies Aq = 27r(5. 

By symmetry, the limit of strong Landau-Zener tran- 




FIG. 2. Landau-Zener tunneling probability P as a func- 
tion of dephasing rate 7 and adiabaticity parameter S. Nu- 
merically exact results (dots) are compared to the prediction 
P — X^n ^2n+i from the jump expansion (black lines), with 
the approximation (19 1 and parameters Aq = 2tt5, Ai — 7r*, 
and r* in (20 1. The deviation does not exceed 0.5% over the 



entire parameter range. 

sitions, Ao(t) 3> Ai(r), yields an inhomogeneous Poisson 
process with rate Ai(t). To obtain Poissonian asymp- 
totics in both Aq and Ai we require an exponential be- 
havior of Pn in Ai as Ai — )■ 00, 



A" 

P — — 



(l + /„(Ai)e-^0 



(18) 



Here, /„(Ai) is a polynomial in Ai, with coefficients de- 
pending on Aq. The first n — 1 of these coefficients are 
determined by the transition probability for 7 = and by 
the asymptotic behavior P„ = 0(A"~^) as Ai — )■ of the 
non-decreasing Markovian (birth) process. By requiring 
that all higher orders of /„ vanish we obtain 



P 



1^ 



r(n,Ai 



{n-l)\ 



{n-l)\ 



1- 



r(n,Ao) 



in 



1)!. 
(19) 

where T[n^x) is the incomplete gamma function. This 
result shows the required asymptotic Poissonian behavior 
of Pn in Aq and Ai. 

Identifying Ai with 7T* , the only free parameter is the 
characteristic Landau-Zener interaction time r*. It can 
be determined by taking the derivative d^Pi at 7 = in 



( 19 ). Pi is the weight of the first order term of the jump 



expansion, for which an exact analytic expression exists 
in terms of the parabolic cylinder functions [29 . With 
the help of their asymptotic expansions |30j we obtain 



7rtanh(5/2(5)y5(l" 
1 - (1 + 27r(5)e- 



-l6<5^ 



-2ttS 



(20) 



Fig.[2]compares the derived approximate transfer prob- 
ability, as follows from ([l9| with r* given by (20 1, to the 



numerically exact value for P((5, 7). One finds a striking 
agreement over the whole parameter range; the devia- 
tions remain below 0.5% and the asymptotic behaviors 



5 



match exactly. This high quahty of the approximation, 
which substantially improves existing ones |19j . makes 
it viable for practical applications such as quantum con- 
trol tasks. Moreover, it suggests that the particular birth 



model (16), (17), which was inspired by the adaptive re- 



summation, captures the essence of the dynamics in this 
problem. 

Conclusion. — The general derivation and the worked 
out examples suggest that the adaptive expansion 
method provides the adequate way of unveiling the in- 
terplay of coherent and incoherent quantum dynamics. 
By guiding our intuition it serves as a natural starting 
point for efficient, analytically tractable approximation 
schemes, which may well advance the purposeful employ- 
ment of incoherent processes. 
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